;***************************************************************
;discription：分析获取植被生长的最适宜温度
;
;param1: nrows(input) 行
;param2: ncols(input) 列
;param3: ndvi_all(input) 一年中12月ndvi
;param4: tmp_all(input)  一年中12月气温
;param5: t_opt(output) 植被生长的最适宜温度
;
;author: wangqiang( wensten )  2015-09-18  chengdu
;***************************************************************
PRO PRO_T_OPTIMAL, nrows, ncols, ndvi_all, tmp_all, t_opt=t_opt
  ;由于模型计算非常复杂，目前取年内最高温度为最适宜温度
  
  ;不知道为啥，第12波段的气温和降水数据有问题
  FOR yindex=0, nrows-1 DO BEGIN
     FOR xindex=0, ncols-1 DO BEGIN 
        ndvi_months = ndvi_all[yindex, xindex, *]
        tmp_months=tmp_all[yindex,xindex, *]
        
        ;当ndvi在一年12个月中达到最高时所在的月份的月气温为植被生长最适宜的温度
        ndvi_max = ndvi_months[yindex, xindex, 0]
        t_opt[yindex, xindex] = tmp_all[yindex, xindex, 0]
        FOR mindex=1, 12-1 DO BEGIN
           ndvi_mon=ndvi_months[yindex, xindex, mindex]
           IF( ndvi_mon GT ndvi_max) THEN BEGIN
              t_opt[yindex, xindex] = tmp_all[yindex, xindex, mindex]
              ndvi_max= ndvi_mon
           ENDIF         
        ENDFOR
     ENDFOR
  ENDFOR
END